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Yale University, New Haven, Connecticut 06511, USA 

We report observations of the b — > d penguin- dominated decays B+ -> K° K+ and B° -> K°K° 
in 316 fb of e + e~ collision data collected with the BABAR detector. We measure the branching 
fractions B(B+ -» F if + ) = (1.61 ± 0.44 ± 0.09) x 10~ 6 and B(B° -> K°K°) = (1.08 ± 0.28 ± 
0.11) x 10" 6 , and the CP- violating charge asymmetry Acp (K°K+) = 0.10 ± 0.26 ± 0.03. Using 
a vertexing technique previously employed in several analyses of all-neutral final states containing 
kaons, we report the first measurement of time-dependent CP- violating asymmetries in B° — > K®K®, 
obtaining S = — 1.28to'.?3 -o/ie and & = ~0.40±0.41±0.06. We also report improved measurements 
of the branching fraction B(B + -» is: 7r + ) = (23.9 ± 1.1 ± 1.0) x 10" 6 and CP-violating charge 
asymmetry Acp (K°tt+) = -0.029 ± 0.039 ± 0.010. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The decays B+ -► K°K+ and B° -> K°K° are ex- 
pected to be dominated by the flavor-changing neutral- 
current process b — > dss, which is highly suppressed in 
the standard model and potentially sensitive to the pres- 
ence of new particles in a way analogous to b — > sss 
decays such as B — > <j>K Q, Assuming top-quark 
dominance in the virtual loop mediating the b — > d tran- 
sition 3], the charge asymmetry in B + — > K°K + and 
the time-dependent CP-violating asymmetry parameters 
in B° — > K®K® are expected to vanish, while contri- 
butions from lighter quarks or supersymmetric particles 
could induce observable asymmetries QJ. It has been 
noted |5J| that the branching fraction and CP asymme- 
tries in B° — ► K°K° are related in a nearly model- 
independent way, providing a sensitive test of the stan- 
dard model description of CP violation. 

In this Letter, we report observations of B + — > K°K + 
and B° — > K°K° using a data sample approximately 
50% larger than the one used in our previous search |fj|. 
(The use of charge-conjugate modes is implied through- 
out this paper unless otherwise stated.) In addition 
to establishing these decay modes, we present measure- 
ments of the time-dependent CP-violating asymmetries 
in B° — > K°K° for the first time. We also report up- 
dated measurements of the branching fraction and charge 
asymmetry in the 5'J7(3)-related decay B + — ► K a ir + . 

The CP asymmetry in B° — > K°K° (observed in the 
K®K® final state) is determined from the difference in 
the time-dependent decay rates for B° and B°. In the 
process e + e~ — > T(AS) — > B°B°, the decay rate /+ (/_) 
is given by Q 

e -|At|/r 

f ± (At) = — [1 ± Ssin(Am rf At) 

4t 

T Ccos(Arn d At)] (1) 

when the second B meson in the event (denoted S tag ) 
is identified as B° (P°). Here At is the time difference 
between the decays of the signal and Ptag mesons, r is 
the average B° lifetime, and Am^ is the B° — B° mixing 
frequency. The amplitude S describes CP violation in the 
interference between mixed and unmixed decays into the 
same final state, while C describes direct CP violation in 
decay. 



The data sample used in this analysis contains (347. 5± 
3.8) x 10 6 T(45) -► BB decays collected by the BABAR 
detector at the Stanford Linear Accelerator Center's 
(SLAC) PEP-II asymmetric-energy e + e~ collider. The 
primary detector elements used in this analysis are a 
charged-particle tracking system consisting of a five-layer 
silicon vertex tracker and a 40-layer drift chamber sur- 
rounded by a 1.5-T solenoidal magnet, and a dedicated 
particle-identification system consisting of a detector of 
internally reflected Cherenkov light. 

We identify two separate event samples corresponding 
to the decays B+ -> K° s h+ and B° -> K°K° S , where 
h + is cither a pion or a kaon. Neutral kaons are recon- 
structed in the mode K® — > 7r + 7r~ by combining pairs of 
oppositely charged tracks originating from a common de- 
cay point and satisfying selection requirements on their 
invariant mass and proper decay time. Candidate h + 
tracks are assigned the pion mass and are required to 
originate from the interaction region and to have a well- 
measured Cherenkov angle (0 C ) consistent with either the 
pion or kaon particle hypothesis. 

For each B° candidate, we require the absolute value 
of the difference AE between its reconstructed energy 
in the center-of-mass (CM) frame and the beam energy 
(y/s/2) to be less than 100 MeV. For B + candidates, we 
require —115 < AE < 75 MeV, where the lower limit 
accounts for an average shift in AE of —45 MeV in the 
K°K + mode due to the assignment of the pion mass to 
the K + candidate. We also define a beam-energy sub- 
stituted mass ibes = \J (s/2 + Pi ■ ps) 2 /P 4 2 — p 2 B , where 
the B-candidate momentum p b and the four- momentum 
of the initial e + e~ state (B;, p,) are calculated in the lab- 
oratory frame. We require 5.20 < ttt-es < 5.29GeV/c 2 for 
B candidates in both samples. To suppress the dominant 
background arising from the process e + e _ — > qq (q = 
u,d,s,c), we calculate the CM angle 9$ between the 
sphericity axis of the B candidate and the spheric- 
ity axis of the remaining charged and neutral particles in 
the event, and require |cos#g| < 0.8. 

After applying all of the above requirements, we find 
2321 (30159) candidates in the B° (B+) sample. The to- 
tal detection efficiencies are given in Tableland include 
the branching fraction for K° s — > 7r + 7r~ [llj and a prob- 
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TABLE I: Summary of results for the total detection efficiencies e, fitted signal yields n, signal-yield significances s (includ- 
ing systematic uncertainty), charge-averaged branching fractions B, and charge asymmetries Acp (including 90% confidence 
intervals). The efficiencies include the branching fraction for — > 7r + 7r~ ' and the probability of 50% for K°K° -> K%K%. 
Branching fractions are calculated assuming equal rates for T{4S) -> B°B° and B + B~ |ljfl. 



Mode 




e (%) 


n 


s(a) 


B (10~ 6 ) 


Acp 


Acp (90%CL) 


B+^ 


K°7T+ 


12.9 ±0.4 


1072 ± 46 ±|? 




23.9 ± 1.1 ± 1.0 


-0.029 ±0.039 ±0.010 


[-0.092, 0.036] 


B+ 


K°K+ 


12.6 ±0.4 


71 ±19 ±4 


5.3 


1.61 ±0.44 ±0.09 


0.10 ±0.26 ±0.03 


[-0.31,0.54] 


B° — > 


K°K° 


8.5 ±0.3 


32 ± 8 ± 3 


7.3 


1.08 ±0.28 ±0.11 







ability of 50% for K°K° -> K°K° 12]. We use data 
and simulated Monte Carlo samples 13] to verify that 
backgrounds from other B decays are negligible. 

A multivariate technique |l4| is employed to determine 
the flavor of the B tag meson in the B sample. Separate 
neural networks are trained to identify primary leptons, 
kaons, low-momentum pions from D* decays, and high- 
momentum charged particles from B decays. Events are 
assigned to one of six mutually exclusive "tagging" cat- 
egories. The quality of tagging is expressed in terms of 
the effective efficiency Q — J2 k £k(l — 2wfc) 2 , where €k 
and Wk are the efficiencies and mistag probabilities, re- 
spectively, for events tagged in category k. We measure 
the tagging performance in a data sample of fully recon- 
structed neutral B decays (Bfl av ) to D^~(-rr + , p + , af), 
where the flavor of the decaying B meson is known, and 
find a total effective efficiency of Q = (30.4 ± 0.3)%. 

The time difference At = Az /(3jc is obtained from the 
known boost of the e + e~ system (f3j = 0.56) and the 
measured distance Az along the beam (z) axis between 
the B° — > K®K® and B tag decay vertices. The posi- 
tion of the Stag vertex is determined from the remaining 
charged particles in the event after removing the four 
tracks composing the signal candidate. Despite the rela- 
tively long lifetime of the K® mesons, the z position of the 
-B-candidate decay point is obtained reliably by exploit- 
ing the precise knowledge of the interaction point using 
the technique described in Ref . [l5j . We compute At and 
its error from a combined fit to the T(4S) — + B°B° decay, 
including the constraint from the known average lifetime 
of the B° meson. Approximately 82% of signal events 
contain a reconstructed from pions that each have 
at least two hits in the silicon vertex tracker, providing 
sufficiently small At uncertainty (0.9 ps) to perform the 
measurement. We require |Ai| < 20 ps and a At < 2.5 ps, 
where a At is the uncertainty on At determined sepa- 
rately for each event. The resolution function for sig- 
nal candidates is a sum of three gaussian distributions 
with parameters determined from the Bfl av sample |ll|. 
The background At distribution has the same functional 
form as the signal resolution function, with parameters 
determined directly from data. 

To obtain the yields and CP violating asymmetry 



parameters in each sample, we apply separate un- 
binned maximum-likelihood fits incorporating discrimi- 
nating variables that account for differences between BB 
and qq events. In addition to the kinematic variables Toes 
and AE, we include a Fisher discriminant T |l6j ] defined 
as an optimized linear combination of the event-shape 
variables J2iP* an d J2iPi c °s 2 6>*, where p* is the CM 
momentum of particle i, 8* is the CM angle between the 
momentum of particle i and the £?-candidate thrust axis, 
and the sum is over all particles in the event excluding 
the B daughters. For the fit to the B + sample we include 
the Cherenkov angle measurement to separate K®tt + and 
K®K + decays. For the B° sample we include At to de- 
termine the CP-violating asymmetry parameters S and 
C simultaneously with the signal yield. 

The likelihood function to be maximized is defined as 
C = exp (- J2i n-i) IIj=i E.^fJi where n r and Vi are 
the yield and probability density function (PDF) for each 
component i in the fit, and N is the total number of 
events in the sample. For the B sample there are two 
components (signal and background) , and the total PDF 
is calculated as the product of the individual PDFs for 
to es , AS, J 7 , and At. The signal At PDF is derived 
from Eq. ^ modified to take into account the mistag 
probability and convolved with the resolution function. 
We combine B + and B~ candidates in a single fit and 
include the PDF for 8 C to determine separate yields and 
charge asymmetries for the two signal components, K®ir 
and K®K, and two corresponding background compo- 
nents. For both signal and background, the K^h^ yields 
are parameterized as n± = n(l =p Acp)/2', we fit directly 
for the total yield n and the charge asymmetry Acp ■ We 
have found correlations among the PDF variables in the 
fit to be negligible in both the B° and B + samples. 

The parameterizations of the PDFs are determined 
from data wherever possible. In both samples, we ex- 
ploit the large sideband regions in toes and AE to de- 
termine all background PDF parameters simultaneously 
with the yields and CP asymmetries in the fits. For the 
B + sample, the large signal Klir + component allows for 
an accurate determination of the peak positions for Toes 
and AE, as well as the parameters describing the shape of 
the PDF for T . The remaining shape parameters describ- 
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ing ttzes and are determined from simulated Monte 
Carlo samples and are fixed in the fit. We use the K^tt + 
parameters to describe signal K Q S K + PDFs in ttieSj AE, 
and T , taking into account the known shift in the mean 
of AE due to the pion-mass hypothesis. For both signal 
and background, the 8 C PDFs are obtained from a sample 
of D* + — > D°ir + (D° — > K~tt + ) decays reconstructed in 
data, as described in Ref. 0|- For the B° sample, all 
shape parameters describing the toes, AE, and T sig- 
nal PDFs are fixed to the values determined from Monte 
Carlo simulation except the peak position for AE, which 
is derived from the results of the fit to the B + sample. 

Several cross-checks were performed to validate the fit- 
ting technique before data in the signal region were ex- 
amined. We checked for biases by performing pseudo- 
experiments where simulated Monte Carlo signal events 
were mixed with background events generated directly 
from the PDFs according to the expected yields in the 
data. The resulting small biases on the yields include ef- 
fects of incorrect particle identification and are accounted 
for in the systematic uncertainties. 

The fit results supersede our previous measurements 
of these quantities and are summarized in Table [I] The 
signal yields for B+ -> K° S K+ and B° -> K° S K° S corre- 
spond to significances of 5.3er and 7.3a (including system- 
atic uncertainties), respectively, and are consistent with 
our previous measurements Q , as well as with recent re- 
sults by the Belle Collaboration [l^. The significances 
are computed by taking the square root of the change in 
21n£ when the appropriate yield is fixed to zero. The 
fit to the B° sample yields S = -1.28±£^ tg'.ie and 
C = — 0.40±0.41±0.06, where the first errors are statisti- 
cal and the second are systematic. The linear correlation 
coefficient between S and C is —32%. 

In Fig. ^ w e compare data and PDFs using the event- 
weighting technique described in Ref. We per- 
form fits excluding the variable being shown; the co- 
variance matrix and remaining PDFs are used to de- 
termine a weight that each event is cither signal (main 
plot) or background (inset). The resulting distributions 
(points with errors) are normalized to the appropriate 
yield and can be directly compared with the PDFs (solid 
curves) used in the fits. We find good agreement be- 
tween data and the assumed shapes in both rriES an d 
AE. In Fig. [2] we display the At distributions for 
KgKg events tagged as B a or B°, and the asymmetry 
A = (N B a — Ngo) I (Nb° + Ng )- The projections are 
enhanced in signal decays by selecting on probability ra- 
tios calculated from the signal and background PDFs (ex- 
cluding At). The likelihood function in the B° -> K°K° 
fit is used to derive Bayesian confidence-level contours 
in the C vs. S plane by fixing (S, C) to specific values, 
refitting the data, and recording the change in — 21og£. 
Figure [3 shows the resulting na contours in the physical 
region defined by S 2 + C 2 < 1 . 

Systematic uncertainties on the signal yields are pri- 
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FIG. 1: Distributions of (left) m E s and (right) AE for sig- 
nal (main plot) and background (inset) (a),(b) 7fs7r + ', (c),(d) 
K®K + , and (e),(f) K®K® candidates (points with error bars) 
in data obtained with the weighting technique described in 
the text. The solid curves represent the assumed shapes used 
in the fits. 
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FIG. 2: Left: 



distributions of At for B° 
50 



K° S K° S decays in 

data tagged as B° (top) or B° (middle), and the asymme- 
try (bottom). The data is enhanced in signal decays using 
requirements on probability ratios. The solid curve repre- 
sents the PDF projection for the sum of signal and back- 
ground, while the dotted curve shows the contribution from 
background only. Right: Likelihood contours in the S vs. C 
plane, where na corresponds to a change in —2 log C of 2.3 for 
n = 1, 6.2 for n = 2, and 11.8 for n = 3. The circle indicates 
the physically allowed region, while the point with error bars 
denotes the result of the fit to data. 



marily due to the imperfect knowledge of the PDF 
shapes. We evaluate this uncertainty by varying the 
PDF parameters that are fixed in the fit within their 
statistical errors, and by substituting different functional 
forms for the PDF shapes. For the charged modes, the 
largest contribution is due to the signal parameterization 
of m ES and AE (3% for K°ir+ , 4% for K° S K+), while 
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for the neutral mode it is due to the potential fit bias 
(8.6%) determined from the pseudo-experiments. We 
use the larger of the value or uncertainty on the back- 
ground asymmetries to set the systematic uncertainty 
on Acp due to potential charge bias |Ti|. We measure 
background asymmetries A C p{K°tt + ) = -0.010 ±0.008 
and Acp{K° s K+) = -0.005 ± 0.009, which arc consis- 
tent with no bias and lead to a systematic uncertainty of 
0.010. The dominant sources of systematic uncertainty 
on S and C are due to the positions of the means in toes 
and AE. The statistical uncertainties of the measured 
values of the CP parameters are in good agreement with 
the expected error values (0.8 ±0.3 for S and 0.6 ±0.2 for 
C), while Monte Carlo studies confirm that the fit tech- 
nique is unbiased for large values of the CP parameters. 

In summary, we have observed the decays B + — > 
K°K + and B° — » K°K° with significances of 5.3<r and 
7.3cr, respectively. The observed branching fractions 
are consistent with recent theoretical estimates 0, Eoj . 
The measured values of the time-dependent CP- violating 
asymmetry parameters in the B° — > K S K® mode re- 
ported here indicate that large positive values of S are 
disfavored, although more data will be needed to confirm 
this result. We have also improved our measurements of 
the branching fraction and CP-violating charge asymme- 
try in B + — > A'!?7r + ; both are consistent with previous 
measurements by other experiments |2l| . 
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